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Abstract. We investigate the operation of pyramidal magneto-optical traps (MOTs) 
microfabricated in silicon. Measurements of the loading and loss rates give insight into 
the role of the nearby surface in the MOT dynamics. Studies of the fluorescence versus 
laser frequency and intensity allow us to develop a simple theory of operation. The 
number of 85 Rb atoms trapped in the pyramid is approximately L 6 , where L < 6 is 
the size in mm. This follows quite naturally from the relation between capture velocity 
and size and differs from the L 3 6 often used to describe larger MOTs. Our results 
constitute substantial progress towards fully integrated atomic physics experiments 
and devices. 
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1. Introduction 



Atom chips are made by micro-fabrication methods and are used to control, trap, and 
manipulate ultra-cold atoms [TJ [2]. Normally, cold atoms are prepared far from the 
surface by a reflection MOT [3j H] or other standard source then transferred to the 
microscopic traps of the atom chip through a delicate sequence of changes to laser beams 
and magnetic fields. These steps make it awkward to load a single trap and difficult 
to load several traps in parallel. We have recently demonstrated a simple alternative 
whereby small pyramids etched into a silicon wafer may be used together with integrated 
current-carrying wires to trap small atom clouds directly from a thermal vapour [5j El [7] . 
Compared with other atom-chip sources this configuration is very simple to operate, and 
because they are micro-fabricated it is easy to form several pyramid MOTs on the same 
chip. 

In this paper we study the operation of integrated MOTs with 85 Rb atoms. In 
Sec. El we describe how the small trapped atom clouds are imaged, then in Sec. 13.11 
we use these images to measure the loading and loss rates of the trap and hence to 
determine the capture velocity. Section [3T21 discusses how the number of atoms is affected 
by the position of the cloud within the pyramid, leading to an understanding of how 
the walls influence the capture of atoms from the vapour and the loss of atoms that are 
already trapped. We note in Sec. 13 .31 that the position of the cloud depends on the laser 
intensity and we measure that effect. Drawing on these results, we are able in Sec. !3.4l to 
understand how the MOT fluorescence signal optimises with respect to laser detuning 
and intensity. Section 13.51 demonstrates that the number of atoms in the cloud scales in 
a simple way as the size of the pyramid is changed. Finally, in Sec. H]we summarise our 
results and briefly discuss future prospects for integrated MOTs. 



2. Experimental Setup 



Figure 1(a) shows a 3mm-thick, 4"-diameter silicon wafer, cut along the [100] crystal 
plane, which we have etched in potassium hydroxide to form concave pyramidal hollows. 
The sides of the pyramids range in length L from 4.2 mm down to 1.3 mm, while the 
apex angle 9 = 2 arctan(l/y / 2) = 70.5° is fixed by the crystal structure. Each set 
of pyramids released from the wafer is coated with 25 nm of aluminium to achieve a 
reflectivity of 75%, so that the pyramid MOTs can work without the need to mask 
any part of the surface [5J. Initial tests were performed on an isolated pyramid, with 
hand-wound coils generating the trapping magnetic field. This was later replaced by 



the chip package shown in Figure 1 (b) containing a 3 x 3 array of pyramids. A carrier 
machined from PEEK supports both the pyramid array and the current- carrying copper 
structure that generates the required array of magnetic quadrupoles. Typical quadrupole 
field gradients used for this package are 0.3 Tin" 1 , produced by a current of 3 A. The 
assembly is mounted in an ultra high vacuum chamber with the pyramid openings 
facing downwards. Rubidium atoms are provided by a current-activated dispenser. 
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(a) The wafer (b) The chip package 



Figure 1. (a) Pyramids etched into a silicon wafer. The wafer is 3mm thick and 4" in 
diameter, (b) A 3 x 3 array of silicon pyramids mounted in PEEK holder 25 x 30 mm 2 . 
The magnetic field required for trapping atoms is provided by a zig-zag array of copper 
wires above and below the pyramids. 



More details relating to the fabrication of these pyramids as well as the chip package 
and laser system are given in previous publications [3 EJ [7]. 

Atoms trapped in the pyramids are detected by imaging their fluorescence onto a 
CCD camera. We use a combination of background subtraction and image filtering to 
distinguish MOT fluorescence from the light scattered by the pyramid surfaces and by 
untrapped rubidium vapour. The MOT is turned on and off by a uniform magnetic field 
of 100 mG along the pyramid axis. This moves the magnetic quadrupole to a position 
where a MOT is not supported, while making almost no change to the fluorescence from 
untrapped atoms. Since the noise in the background-subtracted images is dominated 
by photon shot noise, we are able to improve the visibility of the MOT by averaging 
many images. Figure 2 shows a cloud of 7 x 10 3 atoms. Although we have averaged over 
500 exposures of 40 ms each, the MOT is difficult to distinguish from the background 



in Figure 2(a) The visibility is much improved in Figure 2(b) by convolving Figure 2(a^ 



with a Gaussian [8] having a standard deviation of 5 pixels (45 //m), which is a little 




(a) Before image filtering (b) After image filtering 

Figure 2. A cloud of 7 x 10 3 atoms with Gaussian radius a mot — 100 /xm. a) 
Average of 500 background-subtracted images. The MOT is not clearly visible here 
against the background noise, b) Gaussing filtering of this image clearly reveals the 
MOT. 
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smaller than the ~ 100 /im radius of the MOT. In this way, we can discern MOTs as 
small as 200 atoms at the moment and we anticipate even better sensitivity in future 
pyramids with smoother faces [7j. 



3. Results 



3.1. Loading and loss rates 

The number of atoms in the trap N(t) depends on the balance between the rate R of 
capture from the background Rb vapour and the rate 1/r at which atoms are lost from 
the MOT by collisions with the background gas. The loss due to collisions between 
trapped atoms can be neglected for the traps considered here. Thus the number of 
trapped atoms is well approximated by 

N(t)=N 00 (l-e- t ' T ), (1) 

where Noo = Rt is the steady-state number of atoms. The background Rb pressure 
is adjusted by changing the current in the dispenser. We determine the corresponding 
density and temperature of the Rb vapour, n vap and T vap , by passing a probe laser 
through the chamber and measuring the optical thickness at resonance and the Doppler 
width of the absorption. Figure 3(a) shows N(t) measured at several values of n vap in a 



pyramid of side length L = 4.2 mm. At each density, we turned the trap on at t = and 
recorded 500 MOT images, equally spaced over an interval of 3.3 s. On fitting Eq. (pQ) to 



these measurements, we obtain the loss and capture rates plotted in Figure 3(b) One 
sees that the loss rate is proportional to the Rb vapour density, as expected from the 
relation 1/r = n vap ai oss v vap , where ai oss is the cross section for a trapped atom to be 




0.0 0.2 0.4 0.6 0.8 1.0 1 2 3 4 

t 00 n vap (10 14 m- 3 ) 

(a) Loading curves (b) Loss and capture rates. 



Figure 3. a) Number of atoms in the MOT N(t) versus time for several values 
of the background rubidium vapour density n vap - The lines are least square fits of 
Eq.([T]) to each data set. b) Loss rate 1/r and capture rate R both exhibit a linear 
dependance on rubidium pressure. Linear least-squares fits (solid lines) give 1/r = 
(3.5(2) x lO" 14 ^ + 0.3(2)) s- 1 and R = (9.6(8) x lO" 11 ^ + 1(1)10 3 ) atoms.s- 1 . 
The numbers in parentheses are standard errors given by the fits. 
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kicked out of the trap by collision with a hot Rb atom and v vap is the mean speed of 
atoms in the vapour. At zero background pressure of rubidium, tirj, = 0, our fit to the 
data in Figure 3(b) gives a loss rate of 0.3(2) s -1 , which is consistent with zero. This 
shows that although the walls of the pyramid are close to the trapped atoms, any loss 
to the walls is nevertheless very slow in comparison with the loss rate from collisions 
with the background rubidium vapour. For the measured temperature T vap = 298 K, 
the mean speed is 272 m/s, which yields a measured value for the trap-loss cross section 
of ai oss = 1.29(7) x lCT 16 m 2 . This cross section is related to the velocity needed to 
escape from the trap, given by |9] ai oss = n (4C s /mv esc v vap ) 2 ' 3 , where C 3 = 5.8 x 10~ 48 
Jm 3 is the coefficient of the dominant dipole-dipole interaction and m is the mass of the 
85 Rb atom. Hence the escape velocity is v esc = 2.3(2) ms" 1 . 

The measured capture rate R, shown in Figure 3(b) is also proportional to n vap , 
with a slope of 9.6(8) x 10 _11 m 3 s _1 . In a simple model for this rate, assuming a 
spherical stopping region of cross sectional area A (see Appendix) , this slope is to equal 
8Avc/ (37T 2 v^ ap ), where v c is the maximum velocity of atoms that are captured. We 
can reasonably assume that this capture velocity is twice the escape velocity since the 
distance available for capture is approximately twice the distance needed for escape and 
the friction force is proportional to velocity. Taking v c = 2v esc = 4.6 ms" 1 , we obtain 
a capture area of A = 16 mm 2 , which is reassuringly close to the physical opening area 
of our L = 4.2 mm pyramid. Henceforth we will take it that A ~ I?. Combining 
our results for the capture and loss rates we obtain an expression for the steady state 
number of atoms in the MOT, 

N 00 = Rr^-^— (2) 
3n 2 ai oss v* ap 

Since we expect the capture velocity to be proportional to the size of the MOT, i.e. 
v c oc L, the number of trapped atoms should scale with pyramid size in proportion to 
L 6 , a relation that we verify in Sec. 13.51 



3.2. Position dependent atom number in the MOT 

We are able to adjust the position of the MOT within the pyramid by applying a 
magnetic shim field. FigureH^a) shows a map of the MOT fluorescence versus the shim 
field over a 2.6 mm square region in the horizontal plane 1.8 mm from the apex. When 
the MOT approaches the wall of the pyramid we find that the fluorescence decreases 
quite abruptly due to a decrease of N^. 

In order to determine whether this is due to a decrease in r or reduction of R we set 
up an auxiliary experiment using a larger, aluminium-coated glass pyramid of the same 
apex angle in which we could reliably determine the distance between the MOT and a 
pyramid face by measuring the distance between the MOT and its reflection. The rms 
radius of this cloud was 120(10) /im. Additional loading-curve measurements allowed 
us to determine the loading and loss rates separately, just as in Section 13.11 Over the 
region where the MOT disappears, the loading rate is constant, whereas the loss rate 
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B x (G) d (mm) 

(a) Variation of MOT fluorescence (b) Loss rate to the surface 

Figure 4. a) Variation of MOT fluorescence as the trap is translated horizontally 
within the pyramid by changing a bias magnetic field. The fluorescence decrease when 
the cloud approaches the walls b) Change of loss rate 1/r as a MOT of rms radius 
120 /im approaches a surface. Data: measured rates derived from loading curves. Line: 
Fit to the function Ae~ d ^ 2a ) + 1/tq. Fit parameter a = 119 /im agrees with the size 
of the cloud. Over the same range of distances there is no discernable variation in the 
capture rate R. Inset: Semi- log plot of the same data. 



exhibits a dramatic increase. The loss rate is plotted in Figure@|b) versus the distance 
d from the centre of the MOT to the surface. The data show loss to the surface in 
addition to a background collisional rate 1/t . A good fit is obtained with the function 
1/r = Ae~ d2 ^ 2a2 ^ + 1/tq. This form is motivated by noting that the flux of MOT atoms 
crossing a plane at distance d from the centre is proportional to e~ d2 ^ 2a2 \ when the 
MOT density is a Gaussian of rms radius a. The fit in FigureSJ^b) yields the value 
a = 119/im, in agreement with the measured cloud size. We conclude that atoms are 
lost when the MOT approaches a surface because those at the edge of the cloud collide 
with it. This attenuates the MOT when the cloud comes within within a few hundred 
of the surface. 



Distance from the apex (mm) 
1.50 1.56 1.64 1.71 1.78 1.86 1.93 2.00 2.07 2.14 





1.0 1.2 1.4 1.6 1.8 2.2 2.2 2.4 



Distance from apex (mm) 

Figure 5. (a) Images of the MOT at various distances from the pyramid apex. The 
cloud height is controlled by moving the centre of the magnetic quadrupole using a 
vertical magnetic field. The images are compressed horizontally for the purpose of this 
display, b) Data points: MOT fluorescence integrated over the image and normalised 
to the largest signal. Line: Lorentzian to guide the eye. 
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The bias magnetic field can also move the MOT up and down along the axis of the 
pyramid. Figure[5]^a) shows a series of MOT images taken at different distances from 
the apex of a L = 4.2 mm pyramid. The integrated MOT signals are plotted versus 
distance from the apex in Figure[5^b). We see that the MOT is brightest at distances 
in the range 1.6 — 2.0 mm, centred roughly on 0.43L. Further away from the apex, the 
capture rate R decreases because the MOT beams are cut off by the aperture of the 
pyramid, the last ray (called type-1 in (5J) crossing the axis at 3L/^/32 = 2.2 mm from 
the apex. Beyond this distance a MOT is unable to form. The number of atoms also 
decreases as the MOT moves closer to the apex. This decrease turns on too soon to be 
explained by direct loss of atoms to the surfaces of the pyramid. Instead, we think it 
may be due to a drop in the capture velocity, which causes a decrease in the capture 
rate R. In numerical simulations, we see that atoms being captured need to explore a 
volume as they are slowed down and the available volume is restricted near the apex 
of the pyramid (compared with a plane surface at the same distance), thereby reducing 
the capture velocity. The over-all behaviour of the points in Figure[5](b) is conveniently 
summarised by a Lorentzian function fitted to the data. 

3.3. Intensity- dependent equilibrium position of the MOT 

In a standard 6-beam MOT, pairs of light beams propagate in opposite directions with 
equal intensity and opposite circular polarisation. This produces balanced forces that 
place the MOT at the centre of the magnetic quadrupole, where the magnetic field 
is zero, regardless of the intensity of the light. By contrast, the radiation pressure in 
our 70.5° pyramid is unbalanced vertically at zero magnetic field, due to mismatches 
in both intensity and polarisation of the beams. We therefore expect the equilibrium 
position of the MOT to move up and down on the pyramid axis when the intensity of 
the light is varied. This behaviour is demonstrated by the data in FigureE^a). As the 
intensity increases, the MOT moves further from the apex, indicating that the outward- 




I (mWcm 2 ) 



Figure 6. Distance of the MOT from the apex of the pyramid versus intensity of 
the incident laser beam. As the intensity increases, the MOT moves away from the 
apex because the outward-going beam is more intense. Inset: The movement is small 
compared with the distance to the apex. 
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going light is more intense than the incident beam, despite the fact that it has been 
reflected twice. This is due to a concentration of the reflected light into a smaller cross 
sectional noted in [5]. Over a wide range of intensities, above and below the 

intensity that optimises the MOT, we see a movement that is linear in intensity with 
slope 24.5 /zm/(mW/cm 2 ). As shown by the inset in FigureEl this movement is small 
compared with the total distance to the apex. 

3.4- Understanding the dependence on detuning and intensity 

Figure [3(a) shows the integrated count over the image of the MOT in an L = 4.2 mm 
pyramid versus the detuning 5 of the laser from resonance on the F = 3 — > 4, D2 
transition of 85 Rb at a wavelength of A = 780 nm. The intensity I of the incident 
laser beam is 5 mW / cm 2 , while the repump intensity is 1 mW / cm 2 (the atom number is 
largely insensitive to repump intensity once it exceeds 0.5mW/cm 2 ). We see a peak in 
the fluorescence, centred at 5/2n = —6 MHz and having a width (FWHM) of ~ 5.5 MHz. 
This peak represents the optimum conditions for the MOT and corresponds to having 
7 x 10 3 atoms in the cloud. 

In order to model this behaviour, we consider that the signal should be proportional 
to the number of atoms in the MOT, N^, and to the scattering rate per atom 
S oc 1/(1 + (<5/7) 2 + £1). Here 7 = 3 MHz is half the natural linewidth and £7 
is the saturation parameter that accounts for power broadening. The number is 
proportional to v\ (see Eq. ([2])) and hence to a 4 , where a is the coefficient of friction. 
For this scattering rate S, the corresponding Doppler-cooling friction coefficient is 
a oc 15 /(l + (<5/7) 2 + 6-0 2 - Thus we compare our data with a function / of the form 
N^S oc / = I 5 5 4 / (1 + (S/'j) 2 + £/) 9 . The data points in Fig. [7(a) are well described by 
this function with the choice £ = 2.4, together with a suitable over-all vertical scaling. 




6/2n (MHz) I (mWcm^ 2 ) 



Figure 7. Optimising the MOT fluorescence. Each data point represents the average 
of several experiments, (a) Variation with detuning for an incident intensity of 
5mW/cm . The line is the normalised function / described in the text, (b) Variation 
with incident laser intensity at a detuning of —10 MHz. The dashed (red) line is the 
function /. After correcting for the effect of cloud movement on the fluorescence 
intensity we obtain the solid (blue) curve. 
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The expected value of the coefficient £ can be estimated by considering the 
intensity of the light at the position of the MOT together with its polarisation and 
the corresponding saturation intensity. The total intensity is roughly 61 and the 
polarisation is somewhere between pure circular (I sa t = 1.7mW/cm 2 ) and isotropic 
(I sat = 3.9mW/cni ). Thus, we can reasonably expect £ to be in the range (1.5 — 
3.5) mW'cm 2 , which is entirely consistent with the value that we derive from Fig. |7(a). 
Since we have used the Doppler friction here, the agreement suggests that Doppler 
cooling is an adequate description for this MOT and that the Sisyphus mechanism is 
not playing an important role. One possible reason for the absence of appreciable sub- 
Doppler cooling is the displacement of the cloud from zero magnetic field, which is 
discussed in the previous section. 

Figure [3(b) shows the fluorescence of the same MOT, measured versus incident 
intensity I, with a detuning of S/27T = —10 MHz. The dashed (red) line shows the 
same function / with the same £ = 2.4 and no free parameters, except for the over-all 
scaling. Although this curve captures the general trend, it overestimates the signal at 
high intensity. That is because the MOT fluorescence also changes with the position of 
the cloud (Sec. 13. 2f) and the position changes with the incident intensity (Sec. 13. 3j) . On 
multiplying / by the Lorentzian representing the intensity dependence due to movement 
of the cloud, we obtain the solid (blue) line in Fig. E^b), which is in good agreement 
with the data. 

3. 5. Scaling of atom number with pyramid size 

Throughout this paper we have made use of the scaling law in Eq. [2] that iVoo oc L 2 v^. 
We have further assumed that the MOT is operating within the linear part of the friction 
curve, shown inset in Figure[8l where the frictional force on an atom is proportional to 
its velocity with coefficient a. This gives a capture velocity proportonal to the pyramid 
size L and hence oc L 6 . In order to check this, we measured the MOT fluorescence 
using a range of pyramids having different opening apertures. The smallest pyramids 
were etched in silicon wafers and had openings of L = 4.2,3.5 and 3 mm. We did not 
make the measurement in pyramids smaller than this as the MOT was hard to see and 
difficult to optimise. The etching method was not suitable for making larger pyramids 
because of the very thick wafers and unreasonably long etching times required. In order 
to extend the range to larger apertures, we used an L = 16 mm glass pyramid with the 
same 70.5° geometry and Al coating. An adjustable circular aperture set the diameter 
of the beam entering this larger pyramid to achieve a range of effective sizes, where the 
effective value of L is taken to be the square root of the aperture area. 

The results are shown in Fig. [8] over the range L = 3 — 9 mm, for an incident laser 
intensity of 5mW/cm 2 and a detuning of —8 MHz. We see that up to an aperture size 
of ~ 6 mm, the atom number does indeed follow the L 6 scaling law indicated by the 
solid (red) line. For larger pyramids, the number increases more slowly, and by the 
time the aperture is as large as 1 cm, the scaling is closer to the empirical L 3 6 scaling 
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Figure 8. Scaling law. Red and black data points correspond to MOTs formed in 
silicon pyramids and macroscopic glass pyramids respectively. For pyramids below 
~ 6 mm we find a scaling well described by the relation Noo oc L 6 illustrated by the 
solid (blue) line. The dashed (green) line is the empirical scaling law for large MOTs, 
Noo oc L 3 6 . Inset: Slowing force experiences by atoms versus their velocity. The red 
line shows the approximation F = —av. 



recommended in [10]. This can be understood by considering the capture velocity. 
Scaling linearly from our measurement in Sec. 13. 1[ a 6 mm pyramid has a capture velocity 
of ~ 7m/s. But, as shown in the inset in Fig. [8]this corresponds to the velocity at which 
the frictional force reaches its maximum. Roughly speaking the damping force is linear 
for velocities up to SX, where the Doppler shift is equal to the detuning, beyond which 
the force decreases. Hence, the capture velocity in our experiments grows linearly up to 
a pyramid size of ~ 6 mm, above which it grows ever more slowly. 

4. Summary and outlook 

We have investigated the performance of small pyramid MOTs that are etched into a 
silicon wafer. We have shown how the atom cloud can be imaged against the background 
of scattered light and have used these images to measure the loading and loss rates of 
the trap. By moving the atom cloud within the pyramid, we have measured how the 
number of trapped atoms is affected by proximity to the walls and have shown that the 
main effect is the loss of cold atoms from the tail of the Gaussian density distribution. 
We have studied how the MOT fluorescence varies with position of the cloud, intensity 
of the light and detuning, and we find that the behaviour can be understood in detail 
with the help of a simple Doppler cooling model. The scaling of atom number with 
MOT size has also been understood. This constitutes the first quantitative study of 
MOT behaviour in a small confined space. 
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Miniature pyramids show promise as a chip-scale source of cold atoms for integration 
into other experiments because they are simple and robust, requiring only one circularly 
polarised laser beam and essentially no alignment. For future applications, it would 
be helpful to eliminate the imbalance of beams so that the MOT position remains 
fixed at the zero of magnetic field as the laser intensity varies. This should also 
lead to lower MOT temperatures through the operation of the sisyphus mechanism. 
Such a balance might be achieved by judicious choice of wall coating to lower the 
intensity of the reflected beams. Alternatively, it would be promising if 90° pyramids 
could be fabricated, since these are automatically balanced when the walls have high 
reflectivity [TT1H2]. 

Many experiments require low pressures (~ 1CT 11 mbar) to ensure that the ultracold 
atoms are undisturbed by background gas collisions for many seconds. In such cases, 
the pyramid MOT can be incorporated into a double-sided chip, where both the front 
and back are used. A small vacuum chamber partitioned by the chip would have Rb 
pressure on the pyramid side. Atoms would pass from the MOT to the low pressure side 
through a small aperture etched into the apex of the pyramid, much as in an LVIS source 
[T3| [H] . On arrival, the atoms are trapped and controlled by the usual microfabricated 
wires and optical elements [2]. 

The integrated pyramid MOT is also a promising source of small atom clouds, or 
even single atoms, without the need for extremely high magnetic field gradients |15j . 
From our scaling law studies we conclude that a pyramid with a ~ 1 mm aperture would 
trap a single atom. 



5. Appendix 

In this appendix we derive the expression for the capture rate R used in Sec. 13.11 In 
the Rb vapour, the number density of atoms whose velocity lies within solid angle dfl 
is n vap ^. The Maxwell-Boltzmann probability of a speed between v and v + dv is 
f{v)dv = -^^e~ v2 ^ 2 dv, where (3 = ^j2k B T /m is the most probable speed. Hence the 
flux of atoms moving in this speed range within this solid angle is 

rf$ = n va J^vf{v)dv. (3) 

The number of these per second striking an annulus between b and b+db is 27r6 db d$. Let 
us assume a spherical capture volume with cross sectional area A and radius a = ^ Ajix. 
Integrating over all solid angles and over all speeds up to the capture velocity v c {b) and 
over all impact parameters up to a, the capture rate is 

ra rv c (b) 

R= 2irbdb / n vap vf(v)dv . (4) 
Jo Jo 



Ignoring any deflection of the atoms, the distance available for stopping is 2\J a? — b 2 . 
Since the capture velocity is proportional to this distance, v c (b) = v c ^l — (b/a) 2 , where 
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v c is the maximum capture velocity. After performing the integral in Eq. (J3J) and making 
the approximation [v c j ' (3) 2 <C 1, we obtain the capture rate 

D _ 8Av* 

where the mean speed v vap is equal to ^=/3. 
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